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Abstract. Denote by Q(y— m), with m a square- free positive integer, an imaginary 
quadratic number field, and by O-m its ring of integers. The Bianchi groups are the 
groups SL2(0-m). In the Uterature, there has been so far no example of p-torsion in 
the integral homology of the full Bianchi groups, for p a prime greater than the order of 
elements of finite order in the Bianchi group, which is at most 6. 

However, extending the scope of the computations, we can observe examples of torsion 
in the integral homology of the quotient space, at prime numbers as high as for instance 
p = 80737 at the discriminant —1747. 



1. Introduction 

The Bianchi groups T := SL2(C'_m) may be considered as a key to the study of a larger 
class of groups, the Kleinian groups, which date back to work of Henri Poincare |Poi83| . In 
fact, each non-co-compact arithmetic Kleinian group is commensurable with some Bianchi 
group ||MR03j . A wealth of information on the Bianchi groups can be found in the mono- 
graphs |Fin89| . |EGM98| . |MR03j . In the literature, there has been so far no example of 
p-torsion in the integral homology of the full Bianchi groups, for p a prime greater than 
the order of elements of finite order in the Bianchi group (a recent survey of relevant cal- 
culations has been given in t§enl2| ). In fact, the numerical studies that have been made so 
far, were carried out in the range where the quotient space of hyperbolic 3-space Ti by the 
Bianchi group is often homotopy equivalent to a wedge sum of 2-spheres, 2-tori and Mobius 
bands |,Vog85| . 

We make use of Serre's decomposition |Ser70| of the homology group Hi(r\^; Z,) into the 
direct sum of the free Abelian group with one generator for each element of the class group 
of O-m and the group H™^^(r\^; computed in figures [1] and [2l The first compuations of 
Hi(r; Z) D Hi(r\^; Z) by Swan |Swa71j were on a range of Bianchi groups with vanishing 
cusp-complementary homology H™^''(r\?^; Z). The first example where H™'^''(r\^; is 
non-zero, occurred in an unpublished calculation of Mennicke. Swan's manual computations 
of group presentations have been extended on the computer by Riley [Ril83j : and later 
Vogtmann |Vog85 and Scheutzow |Sch92| systematically computed H^"**''(r\'H; Q) for a 



large range of Bianchi groups. But they were still in in the range where r\^ admits no 
homological torsion. Aranes [Aral O] has computed cell complexes for the Bianchi groups for 
all m ^ 100, and Yasaki |YaslO) has obtained GL2(C'_m)-cell complexes (with the Voronoi 
model) for the same range as well as all cases where O is of class number 1 or 2. This includes 
two cases, m = 74 and m = 86, where some 2-torsion appears in H™*^^(r\'H; Z), but the 
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latter two authors have not yet provided homology computations. When the absolute value 
of the discriminant gets greater, torsion in the integral homology of the quotient space 
appears (see figure [2]) at prime numbers as high as for instance 80737 at the discriminant 
— 1747, whereas the order of elements of finite order in SL2(0-m) is at most 6. A growth of 
the torsion in the Abelianization of the Bianchi groups with respect to the covolume can be 
observed, which is in concordance with the predictions of |BV12| . We can also observe that 
the occurring torsion subgroups are quite likely to occur as squares, but this is no general 
principle, because the discriminant —431 produces a counterexample to this phenomenon. 

In order to obtain the results of figures [T] and [21 in section [7] we fill out Swan's concept 
|Swa7 1] and elaborate algorithms to compute a fundamental polyhedron for the action of 
the Bianchi groups on hyperbolic 3-space. Other algorithms based on the same concept 
have independently been implemented by Cremona |Cre84| for the five cases where O-m 
is Euclidean, and by his students Whitley |Whi90) for the non-Euclidean principal ideal 
domain cases, Bygott |Byg98| for a case of class number 2 and Lingham ( |Lin05| . used 
in |CL07| ) for some cases of class number 3; and finally Aranes |AralO| for arbitrary class 
numbers. The algorithms presented in subsection [7] come with an implementation |RahlO] 
for all Bianchi groups; and we make explicit use of the cell complexes it produces. The 
provided implementation |RahlO| has been validated by the project PLUME of the CNRS, 
and is subject to the certificate C3I of the GENCI and the CPU. Other results obtained 
with the employed implementation are described in [Rahllj and |RS12| . On the computing 
clusters of the Weizmann Institute of Science, this implementation has been applied to 
establish a database of cell complexes for over 180 Bianchi groups, using over fifty processor- 
months. This database includes all the cases of ideal class numbers 3 and 5, most of the 
cases of ideal class number 4 and all of the cases of discriminant absolute value bounded 
by the number 500. 

A computational advantage is the shortcut that we obtain in section H] by linking the 
Borel-Serre compactification of the quotient space with Floge's compactification in a long 
exact sequence, based on the recent paper [Rahl2b| . Floge's compactification admits a 
computationally easier cell structure, and we can explicitly calculate the equivariant Leray- 
Serre spectral sequence associated to it. In section O we describe how to assemble the 
homology of the Borel-Serre compactified quotient space and the Farrell cohomology of a 
Bianchi group to its full group homology with trivial Z-coefficients. Here, we divide by the 
center of SL2(C'_m), consisting of plus and minus the identity matrix, yielding PSL2(C'_m)- 
As the center of SL2(C'_m) is the kernel of its action on hyperbolic 3-space, this does not 
change the quotient space. And for T := PSL2(C'_m), general formulae for its Farrell 
cohomology have been given [Rahl2| (based on [Rahllbj ). 

1.1. Organization of the paper. We print the isomorphism types of H™'^^(r\^; Z) that 
were obtained in figures [1] and [2l The homology group Hi(r\^; is a direct sum of the 
former and the free Abelian group with rank the cardinality of the class group of O-m, 
which we also print. There is an inclusion of Hi(r\'H; Z) into the group homology Hi(r; Z); 
and the latter group homology is a quotient of the direct sum of Hi(r\^; Z) and the Farrell 
supplement that has been computed and printed in a separate column. In section [21 we 
define the Bianchi fundamental polyhedron, which induces our cell structure on r\^- We 
use it in section [3] to obtain the Floge cellular complex, which we connect in section [4] to the 
Borel-Serre compactification of r\^- Then we proceed to Hi(PSL2(C'_m); ^) in section [5l 
describe Swan's concept in section [6] and its realization in section [T] 
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Figure 1. The cusp-complementary homology H^*^*^^ (r\W; ^) for the ab- 
solute values of the discriminant A fulfilling |A| ^ 415. 
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Figure 2. H^"^^(r\^; Z), with its torsion decomposed into prime power 
factors, for some greater absolute values of the discriminant. 
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2. The Bianchi fundamental polyhedron 



Let m be a squarefree positive integer and consider the imaginary quadratic number field 
Qiy^—m) with ring of integers O-m, which we also just denote by O. Consider the familiar 
action by fractional linear transformations (we give an explicit formula for it in lemma [2T]) 
of the group T := SL2(C') C GL2(C) on hyperbolic three-space, for which we will use the 
upper-half space model 7^. As a set, 



The Bianchi-Humbert theory |Bia92| . |Huml5| gives a fundamental domain for this ac- 
tion. We will start by giving a geometric description of it, and the arguments why it is a 
fundamental domain. 

Definition 1. A pair of elements (^, A) G is called unimodular ii the ideal sum fj,0 + XO 
equals O. 

The boundary of Ti is the Riemann sphere dTi = C U {oo} (as a set), which contains the 
complex plane C. The totally geodesic surfaces in Ti are the Euclidean vertical planes (we 
define vertical as orthogonal to the complex plane) and the Euclidean hemispheres centred 
on the complex plane. 

Notation 2. Given a unimodular pair (/i. A) e with /U 7^ 0, let S^^x C Ti denote the 
hemisphere given by the equation \fiz — Ap -|- |/ipC^ = 1- 



This hemisphere has centre A/^ on the complex plane C, and radius Let 
B := {{z,C) e 'H: The inequality \^iz - Ap + \n\'^C^ ^ 1 



Then B is the set of points in Ti which lie above or on all hemispheres S^^\. 

Lemma 3 ( |Swa71| ). The set B contains representatives for all the orbits of points under 
the action o/SL2(C') on %. 

The action extends continuously to the boundary which is a Riemann sphere. 
In r := SL2(C'-m), consider the stabiliser subg roup of the point oo G dl~L~ In the cases 
m = 1 and m = 3, the latter group contains some rotation matrices like i^J-^ ^^o^ ' ^^^^^ 
we want to exclude. These two cases have been treated in |Men79| . |SV83| and others, and 
we assume m 7^ 1, m ^ 3 throughout the remainder of this article. Then, 



which performs translations by the elements of O with respect to the Euclidean geometry 
of the upper-half space %. 

Notation 4. A fundamental domain for Too in the complex plane (as a subset of 8%) is 
given by the rectangle 



'H = {{zX) gCxM|C>0}. 



is fulfilled for all unimodular pairs (/i, A) G with /i / }. 





{x + y^/^ gC|0^x^1, O^ysCl}, 




m = 1 or 2 mod 4, 
m = 3 mod 4. 
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And a fundamental domain for Too in % is given by 

^oo :={(z,C) eT^lzGDo}. 
Definition 5. We define the Bianchi fundamental polyhedron as 

D:=D^n B. 

It is a polyhedron in hyperbolic space up to the missing vertex oo, and up to missing 
vertices at the singular points if O is not a principal ideal domain (see subsection 16. 2p . As 
Lemma [3] states T ■ B = %, and as Too • = H yields Too ■ D = B, we have T ■ D = T-L. 
We observe the following notion of strictness of the fundamental domain: the interior of 
the Bianchi fundamental polyhedron contains no two points which are identified by F. 
Swan proves the following theorem, which implies that the boundary of the Bianchi funda- 
mental polyhedron consists of finitely many cells. 

Theorem 6 ([Swa71]). There is only a finite number of unimodular pairs (A, /u) such that 
the intersection of S^^x with the Bianchi fundamental polyhedron is non-empty. 

Swan further proves a corollary, from which it can be deduced that the action ofT onTi 
is properly discontinuous. 

3. The Floge cellular complex 

In order to obtain a cell complex with compact quotient space, we proceed in the 
following way due to Floge |Flo83| . The boundary of Ti is the Riemann sphere dTi, 
which, as a topological space, is made up of the complex plane C compactified with 
the cusp oo. The totally geodesic surfaces in Ti are the Euclidean vertical planes (we 
define vertical as orthogonal to the complex plane) and the Euclidean hemispheres centred 
on the complex plane. The action of the Bianchi groups extends continuously to the 
boundary dTi. Consider the cellular structure on Ti induced by the F-images of the Bi- 
anchi fundamental polyhedron. The cellular closure of this cell complex in TiL) 8% con- 
sists of n and (Q(^/^) U {oo}) C (C U {oo}) = dU. The SL2(C'_m)-orbit of a cusp ^ in 
(Q(\/— m) U {oo}) corresponds to the ideal class [(A,/x)] of O-m- It is well-known that this 
does not depend on the choice of the representative ^. We extend our cell complex to a 

cell complex X by joining to it, in the case that is not a principal ideal domain, the 

SL2(C'_m)^orbits of the cusps ^ for which the ideal (A,//) is not principal. At these cusps, 

we equip X with the "horoball topology" described in |Flo83| . This simply means that the 
set of cusps, which is discrete in 8%^ is located at the hyperbolic extremities of X : No 
neighbourhood of a cusp, except the whole X, contains any other cusp. 

We retract X in the following, SL2( O-m )~equi variant, way. On the Bianchi fundamental 
polyhedron, the retraction is given by the vertical projection (away from the cusp oo) onto 
its facets which are closed in 7^ U 8%. The latter are the facets which do not touch the 
cusp oo, and are the bottom facets with respect to our vertical direction. The retraction 
is continued on % by the group action. It is proven in |Flo80| that this retraction is 
continuous. We call the retract of X the Floge cellular complex and denote it by X. So 
in the principal ideal domain cases, X is a retract of the original cellular structure on %, 
obtained by contracting the Bianchi fundamental polyhedron onto its cells which do not 
touch the boundary oi %. In [RFllj . it is checked that the Floge cellular complex is 
contractible. 
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4. Connecting Floge cell complex and Borel-Serre compactification 

Let r be a Bianchi group with O admitting as only units {±1}, i.e. we suppose O not to 
be the Gaussian or Eisenstein integers. In the latter two cases, the problem of the singular 
cusps treated here does not occur anyway. Let Tj be the torus attached at the cusp i of F, 
and let Xi and yi denote the cycles generating Hi(Tj). Let P be the Bianchi fundamental 
polyhedron of F. Write "hyp. cells" for cells in the interior of hyperbolic space. Consider 
the short exact sequence of chain complexes obtained from collapsing the singular tori, 



singular 

0— e (T, 







singular 



— ^ {xs,ys 



{P}- 

\d3 



\d3 



a any cusp 

(T,) e (hyp. 2-cells) 



^ any cusp 



0^2 



(Too) e (hyp. 2-cells) 



{xc, Vc) © (hyp. edges) (xoo, yoo) © (hyp. edges) ► 

e ai e <9i 



any cusp any cusp 

" ► (c) © (hyp. vertices) ► (c) © (hyp. vertices) — " 0. 

c c 



Poincare's theorem on fundamental polyhedra tells us that ds^P) 



any cusp 

U T, 



and 



hence (?3(-P) = (Too)- From |Rahl2b| . we see that for every cusp c, there is a chain of 
hyperbolic 2-cells that we denote by ch{xc) and which is mapped to the cycle Xc by 82- 
And furthermore, yc is in the cokernel of 82 (of course, this holds up to the appropriate 
permutation of the labels Xc and yc)- This implies that d2{ch(xoo)) = x^o and y^o is in 
the cokernel of 82- As the quotient space is path- wise connected, the cokernel of 81 is 
isomorphic to Z. The above information tells us that the long exact sequence induced on 
integral homology by the map /? concentrates in 



singular 

(T. 



/32 



any cusp 

( (t,))/<u.t.>©h: 



cusp 
2 



-*Hrp0,(cMx,)) 




where the maps without labels are the obvious restriction maps making the sequence 
exact; and where H™^^ and H™^'' are generated by cycles from the interior of r\^- 



any cusp 



Note that H™*^^ 0^(c/i(xs)) is non-naturally isomorphic to {'^c) j / (Ucfc)®^™^^ ^ 

namely collapsing a torus moves its 2-cycle into a bubble ch{xs) emerging adjacent to 
the singular cusp s in the Floge complex. 
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5. The equivariant spectral sequence to group homology 

Let r := PSL2(C'_m), and let X be the Floge cellular complex of section [3l the cell 
structure of which we subdivide until the cells are fixed pointwise by their stabilisers. We 
describe now how to assemble the homology of the Borel-Serre compactified quotient space 
(issue of the previous section) and the Farrell cohomology of F, for which general formulae 
have been given in |R,ahl2| (based on |R,ahllb) V to the full group homology of F with trivial 
Z-coefficients. We proceed following |Bro82| VII] and |SV83| . Let us consider the homology 
H^,(F;C,(X)) of F with coefficients in the cellular chain complex C,{X) associated to X] 
and call it the F-equivariant homology of X. As X is contractible, the map X — )• pt. to the 
point pt. induces an isomorphism 

H,(F; C,{X)) ^ H,(F; C,{pt.)) ^ H,(F; Z). 

Denote by X^ the set of p-cells of X, and make use of that the stabiliser Fq- in F of any 
p-cell (T oi X fixes a pointwise. Then from 

Cp{X)= Indf^Z, 

Shapiro's lemma yields 

H,(F;Cp(X))^ H,(F,;Z); 

and the equivariant Leray/Serre spectral sequence takes the form 

4,7= H,(F.; Z) ^ Hp+,(F; C.(X)), 

converging to the F-equivariant homology of X, which is, as we have already seen, iso- 
morphic to Hp+q(F; Z) with the trivial action on the coefficients Z. 

As in degrees above the virtual cohomological dimension, which is 2 for the Bianchi 
groups, the group homology is isomorphic to the Farrell cohomology, we obtain the iso- 
morphism type from the above mentioned general formulae. 

In the lower degrees q £ {0,1,2}, the following terms remain on the ii^^-page, which is 
concentrated in the columns p = 0, 1, 2: 

q = 2 © 2i © 2-torsion © 3-torsion 2-torsion © 3-torsion 

s singular 



q = 1 © ffi Farrell supplement 2-torsion © 3-torsion 

s singular 

q = Z lii{r\X;Z) R-jirXX; Z) 

where the "Farrell supplement" is the cokernel of the map 

Hi(F„;Z)3^ Hi(F„;Z). 

induced by inclusion of finite cell stabilisers. As the cells are fixed pointwise by their 
stabilisers, we see that for q > 0, the E'p ^^-terms are concentrated in the two columns p = 
and p = 1. We compute the bottom row {q = 0) of the above spectral sequence as the 
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homology of the quotient space r\^- Then we infer from section |4] that the rational rank 
of the differential q is the number of non-trivial ideal classes of O—'^yi. 

Using Serre's decomposition of the homology group Hi(r\^; into the direct sum of 
the free Abelian group with one generator for each element of the class group of O-m 
and the group H^'^''^(r\'H; Z), and using the long exact sequence of section HI we see that 
Hi(r\X; Z) ^ Rf^{r\n; Z)®Z. This has made it possible ot compute Rf^{r\n; Z) 
from the quotient space of the Floge cellular complex in figures [T] and [2J Furthermore, we 
have an inclusion of Hi(r\'H; Z) into the group homology Hi(r; Z); and the latter group 
homology is a quotient of the direct sum of Hi(r\'H; Z) and the Farrell supplement. 



6. Swan's concept to determine the Bianchi fundamental polyhedron 

This section recalls Richard G. Swan's work |Swa71| . which gives a concept — from the 
theoretical viewpoint — for an algorithm to compute the Bianchi fundamental polyhedron. 
The set B which determines the Bianchi fundamental polyhedron has been defined using 
infinitely many hemispheres. But we will see that only a finite number of them are sig- 
nificant for this purpose and need to be computed. We will state a criterion for what is 
an appropriate choice that gives us precisely the set B. This criterion is easy to verify in 
practice. Suppose we have made a finite selection of n hemispheres. The index i running 
from 1 through n, we denote the i-th hemisphere by S{ai), where is its centre and given 
by a fraction Oj = ^ in the number field Q{-^—m). Here, we require the ideal {Xi,^i) to 
be the whole ring of integers O. This requirement is just the one already made for all the 
hemispheres in the definition of B. Now, we can do an approximation of notation [21 using, 
modulo the translation group Tqq, a finite number of hemispheres. 

Notation 7. Let i3(ai, . . . , a„) := {{z,C) G 'H: The inequality l/xz - Ap + |/ipC^ ^ 1 
is fulfilled for all unimodular pairs (/x. A) G with ^ = Oi -|- 7, for some i G {1, . . . ,n} 

and some 7^0}. Then B{ai, . . . , an) is the set of all points in T-L lying above or on all 
hemispheres S{ai -|- 7), z = 1, . . . , n; for any 7 G O. 

The intersection -B(ai, . . . , a„)nDoo with the fundamental domain for the translation 
group Too, is our candidate to equal the Bianchi fundamental polyhedron. 

6.1. Convergence of the approximation. We will give a method to decide when 

B[ai, . . . ,an) = B. This gives us an effective way to find B by adding more and more 
elements to the set {ai, . . . , a„} until we find B{ai, . . . , q„) = B. We consider the boundary 
dB{ai, . . . , On) of B{ai, . . . , a„) mHuC It consists of the points (^;, C) G 'HUC satisfying 
all the non-strict inequalities l/iz— Ap-|-|/ipC^ ^ 1 that we have used to define i?(ai, . . . , a„), 
and satisfy the additional condition that at least one of these non-strict inequalities is an 
equality. We will see below that dB{ai, . . . ,an) carries a natural cell structure. This, 
together with the following definitions, makes it possible to state the criterion which tells 
us when we have found all the hemispheres relevant for the Bianchi fundamental polyhedron. 

Definition 8. We shall say that the hemisphere 5"^^;^ is strictly below the hemisphere S"^ q, 
at a point z G C if the following inequality is satisfied: 

a ^ _ J_ A ^ 1 
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This is, of course, an abuse of language because there may not be any points on Sp^a or 
S^^x with coordinate z. However, if there is a point [z^Q) on S^^\, the right hand side of 
the inequahty is just — C^- Thus the left hand side is negative and so of the form — (C')^- 
Clearly, (z, C,') G Sp^a and C' > We will further say that a point (z, C) G U C is strictly 
below a hemisphere Sfj^^x, if there is a point {z,(') G 'S'^.A with > ^. 

6.2. Singular points. We call cusps the elements of the number field K = Q(-^/^^m ) con- 
sidered as points in the boundary of hyperbolic space, via an embedding K C CU {00} = dTi 
We write 00 = g, which we also consider as a cusp. It is well-known that the set of cusps is 
closed under the action of SL2(C') on dH; and that we have the following bijective corres- 
pondence between the SL2(C')-orbits of cusps and the ideal classes in O. A cusp ^ is in the 

SL2(0)-orbit of the cusp if and only if the ideals (A', /x') and (A, fi) are in the same ideal 
class. It immediately follows that the orbit of the cusp 00 = ^ corresponds to the principal 
ideals. Let us call singular the cusps ^ such that (A, /x) is not principal. And let us call 
singular points the singular cusps which lie in dB. It follows from the characterisation of 
the singular points by Bianchi that they are precisely the points in C C dH which cannot be 
strictly below any hemisphere. In the cases where O is a principal ideal domain, K U {00} 
consists of only one SL2(C)-orbit, so there are no singular points. We use the following 
formulae derived by Swan, to compute representatives modulo the translations by T^o, of 
the singular points. 

Lemma 9 f |Swa71) ). The singular points of K, mod O, are given by 
P^I±^E^^ where p,r, s G Z, s > 0, ^ < r ^ |, s"^ ^ r"^ + m, and 

• if m = 1 or 2 mod 4, 

s ^ 1, s\r'^ + m, the numbers p and s are coprime, and p is taken mod s; 

• if m = 3 mod 4, 

s is even, s ^ 2, 2s \ r'^ + m, the numbers p and | are coprime; p is taken mod |. 

The singular points need not be considered in Swan's termination criterion, because they 
cannot be strictly below any hemisphere ^^^a- 

6.3. Swan's termination criterion. We observe that the set of z G C over which some 
hemisphere is strictly below another is C or an open half-plane. In the latter case, the 
boundary of this is a line. 

Notation 10. Denote by L(|, ^) the set of z G C over which neither Sp^a is strictly below 
S^^x nor vice versa. 

This line is computed by turning the inequality in definition [8] into an equation. Swan 
calls it the line over which the two hemispheres agree, and we will see later that the most 
important edges of the Bianchi fundamental polyhedron lie on the preimages of such lines. 
We now restrict our attention to a set of hemispheres which is finite modulo the translations 
in Too- Consider a set of hemispheres S{ai + 7), where the index i runs from 1 through n, 
and 7 runs through O. We call this set of hemispheres a collection, if every non-singular 
point z £ C C &H is strictly below some hemisphere in our set. Now consider a set 
B{ai, . . . , a„) which is determined by such a collection of hemispheres. 

Theorem 11 (Swan's termination criterion |Swa71| ). We have B{ai, . . . ,an) = B if and 
only if no vertex of dB{ai, . . . , an) can he strictly below any hemisphere Sfj,^x- 
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In other words, no vertex v of dB{ai, . . . , q„) can lie strictly below any hemisphere S^^\. 
Let us call height the coordinate of the upper-half space model introduced at the beginning 
of section [21 With this criterion, it suffices to compute the cell structure of dB{ai, . . . , a„) 
to see if our choice of hemispheres gives us the Bianchi fundamental polyhedron. This has 
only to be done modulo the translations of Too, which preserve the height and hence the 
situations of being strictly below. Thus our computations only need to be carried out on a 
finite set of hemispheres. 

6.4. Computing the cell structure in the complex plane. We will in a first step 
compute the image of the cell structure under the homeomorphism from dB{ai, . . . , a„) to 
C given by the vertical projection. For each 2-cell of this structure, there is an associated 
hemisphere 5"^, a- The interior of this 2-cell consists of the points z G C, where all other 
hemispheres in our collection are strictly below 5*^, a- Swan shows that this is the interior of 
a convex polygon. The edges of these polygons lie on real lines in C specified in notation 1101 

A vertex is an intersection point z of any two of these lines involving the same hemisphere 
S^^x, if a-ll other hemispheres in our collection are strictly below, or agree with, a at z. 

6.5. Lifting the cell structure back to hyperbolic space. Now we can lift the cell 
structure back to dB(ai, . . . , On), using the projection homeomorphism onto C. The preim- 
ages of the convex polygons of the cell structure on C, are totally geodesic hyperbolic 
polygons each lying on one of the hemispheres in our collection. These are the 2-cells of 
dB{ai, . . .,«„)■ 

The edges of these hyperbolic polygons lie on the intersection arcs of pairs of hemispheres 
in our collection. As two Euclidean 2-spheres intersect, if they do so non-trivially, in a circle 
centred on the straight line which connects the two 2-sphere centres, such an intersection arc 
lies on a semicircle centred in the complex plane. The plane which contains this semicircle 
must be orthogonal to the connecting line, hence a vertical plane in Ti. We can alternatively 
conclude the latter facts observing that an edge which two totally geodesic polygons have 
in common must be a geodesic segment. Lifting the vertices becomes now obvious from 
their definition. This enables us to check Swan's termination criterion. 

We will now sketch Swan's proof of this criterion. Let P be one of the convex polygons of 
the cell structure on C. The preimage of P lies on one hemisphere S{ai) of our collection. 
Now the condition stated in theorem [TT] says that at the vertices of P, the hemisphere 
S{ai) cannot be strictly below any other hemisphere. The points where S{ai) can be 
strictly below some hemisphere constitute an open half-plane in C, and hence cannot lie 
in the convex hull of the vertices of P, which is P. Theorem [TT] now follows because C is 
tessellated by these convex polygons. 
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7. Algorithms realizing Swan's concept 

From now on, we will work on putting Swan's concept into practice. We can reduce the 
set of hemispheres on which we carry out our computations, with the help of the following 
notion. 

Definition 12. A hemisphere S^^\ is said to be everywhere below a hemisphere Sp^a when: 



A a 

Note that this is also the case when S^^\ 



1 1 

¥\~\iA- 

= Sjs^a- Any hemisphere which is everywhere 



below another one, does not contribute to the Bianchi fundamental polyhedron, in the 
following sense. 

Proposition 13. Let 5(a„) be a hemisphere everywhere below some other hemisphere 

S{ai), where i G {1, . . . , n — 1}. 

Then B{ai, . . . , q„) = B{ai, . . . , q„_i). 

Proof. Write ctn = ji and Qj = ^ with A, fJ.,0,T G O. We take any point (z, () strictly below 
S^^x and show that it is also strictly below Sr,9- In terms of notation [71 this problem looks 
as follows: we assume that the inequality \fiz — Ap + < 1 is satisfied, and show that 

this implies the inequality \tz — 9\'^ + |rp(^^ < 1. The first inequality can be transformed 
into 

Hence, 



1 



— —\ + O < TT- We will insert this into the triangle 
inequality for the Euclidean distance in C x M applied to the three points {z, (), (^,0) and 
(-,0), which is 



So we obtain 



e|2 



< 



\_e_ 

fl T 

1 





A 




z 







+ 



By definition [121 the expression on the 



right hand side is smaller than or equal to Therefore, we take the square and obtain 

I fl|2 



+ < |7p ) which is equivalent to the claimed inequality. 
Another notion that will be useful for our algorithm, is the following. 



□ 



Definition 14. Let z G C be a point lying within the vertical projection of 5*^, a- Define 
the lift on the hemisphere 5^^^ of ^ as the point on S^^^x the vertical projection of which 
is z. 

Notation 15. Denote by the hemisphere list a list into which we will record a finite number 
of hemisphere s S{ai), . . ., S{an)- Its purpose is to determine a set B{ai, . . . , an) in order 
to approximate, and finally obtain, the Bianchi fundamental polyhedron. 
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7.1. The algorithm computing the Bianchi fundamental polyhedron. We now de- 
scribe the algorithm that we have reahzed using Swan's description; it is decomposed into 
algorithms [T] through |3] below. 

Initial step. We begin with the smallest value which the norm of a non-zero element 
/i G O can take, namely 1. Then is a unit in O, and for any X € O, the pair (/i, A) 
is unimodular. And we can rewrite the fraction 7- such that u = 1. We obtain the unit 
hemispheres (of radius 1), centred at the imaginary quadratic integers X £ O. We record 
into the hemisphere list the ones which touch the Bianchi fundamental polyhedron, i.e. the 
ones the centre of which lies in the fundamental rectangle Dq (of notation for the action 
of Too on the complex plane. 

Step A. Increase \^\ to the next higher value which the norm takes on elements of O. 
Run through all the finitely many /i which have this norm. For each of these /i, run through 
all the finitely many A with ^ in the fundamental rectangle Dq. Check that (//, A) = O and 
that the hemisphere S^^^x is not everywhere below a hemisphere Sj^^a in the hemisphere list. 
If these two checks are passed, record {fi, A) into the hemisphere list. 

We repeat step A until has reached an expected value. Then we check if we have 
found all the hemispheres which touch the Bianchi fundamental polyhedron, as follows. 

Step B. We compute the lines L{^, ^) of definition [TOl over which two hemispheres 
agree, for all pairs 5^^q,, S'^^a hi the hemisphere list which touch one another. 
Then, for each hemisphere Sp^ai we compute the intersection points of each two lines L{^, ^) 
and L(^, ^) referring to ^. 

We drop the intersection points at which 5^^q, is strictly below some hemisphere in the list. 
We erase the hemispheres from our list, for which less than three intersection points remain. 
We can do this because a hemisphere which touches the Bianchi fundamental polyhedron 
only in two vertices shares only an edge with it and no 2-cell. 

Now, the vertices of B(ai, . . . , a„) Pi D^o are the lifts of the remaining intersection points. 
Thus we can check Swan's termination criterion (theorem llip . which we do as follows. We 
pick the lowest value C > ^01 which (z, £ Ti is the lift inside Hyperbolic Space of a 
remaining intersection point z. 

If C ^ then all (infinitely many) remaining hemispheres have radius equal or smaller 
than so (z, Q cannot be strictly below them. So Swan's termination criterion is fulfilled, 
we have found the Bianchi fundamental polyhedron, and can proceed by determining its 
cell structure. 

Else, ^ becomes the new expected value for We repeat step A until reaches ^ and 
then proceed again with step B. 



14 



ALEXANDER D. RAHM 



Algorithm 1 Computation of tlie Bianclii fundamental polyliedron 
Input: A square- free positive integer m. 

Output: The hemisphere list, containing entries ^(q;!),. . . ,S'(q„) such that 
B{ai, ...,an) = B. 



Let O be the ring of integers in Q(a/— m). 

Let ho be the class number of O. Compute ho. 

Estimate the highest value for |^| which will occur in notation [7] by 



the formula E 



^"'/io-2m + i m = 3mod4, 



1 2\mho — 19m, else. 
M := 1. 

Swan's_cancel_criterion_fulfilled := false. 

while Swan's_cancel_criterion_fulfilled = false, do 
vi^hile TV ^ £^ do 

Execute algorithm [2] with argument M. 

Increase M to the next greater value in 

the set {y^nFm + p \ G N} of values of the norm on O. 
end while 

Compute with algorithm [3j 
if jj- , then 

All (infinitely many) remaining hemispheres have radius 
smaller than , 

so {z,C) cannot be strictly below any of them. 
Swan's_cancel_criterion_fulfilled :— true, 
else 

^ becomes the new expected lowest value for jj- : 

end if 
end while 



Proposition 16. The hemisphere list, as computed by algorithm^ determines the Bianchi 
fundamental polyhedron. This algorithm terminates within finite time. 

Proof. 

• The value is the minimal height of the non-singular vertices of the cell complex 
dB{ai, . . . , On) determined by the hemisphere list {S{ai), . . . , S{an)}. 

All the hemispheres which are not in the list, have radius smaller than jj-. 
By remark [181 the inequality C ^ will become satisfied; and then no non- 
singular vertex of dB{ai, . . . ,a„) can be strictly below any of them. Hence by 
theorem [m B{ai, . . . ,a„) = B; and we obtain the Bianchi fundamental polyhed- 
ron as B{ai, . . . , On) n Dqo- 

• We now consider the run-time. By theorem [H the set of hemispheres 

{5"^, A I 5"^, A touches the Bianchi Fundamental Polyhedron} 
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is finite. So, there exists an Sf^^x for which the norm of fi takes its maximum on 
this finite set. The variable J\f reaches this maximum for |^| after a finite number 
of steps; and then Swan's termination criterion is fulfilled. The latter steps require 
a finite run-time because of propositions [19] and [20l 

□ 

Swan explains furthermore how to obtain an a priori bound for the norm of the fi G O 
occurring for such hemispheres S^^a- But he states that this upper bound for |^| is much 
too large. So instead of the theory behind theorem [6l we use Swan's termination criterion 
(theorem [TT] above) to limit the number of steps in our computations. We then get the 
following. 

Observation 17. We can give bounds for in the cases where K is of class number 1 
or 2 (there are nine cases of class number 1 and eighteen cases of class number 2, and we 
have done the computation for all of them). They are the following: 



|A|+1 



K of class number 1: l/^l ^ 2 
f|/i| ^ 3|A|, 

K of class number 2: ^ (5 ^ ^)\A\, else 



where A is the discriminant of = Q(\/— m), i.e., |A| 




4m, else. 



Remark 18. In algorithm [H we have chosen the value E by an extrapolation formula for 
observation 1171 If this is greater than the exact bound for the algorithm computes 
additional hemispheres which do not contribute to the Bianchi fundamental polyhedron. 
On the other hand, if E is smaller than the exact bound for |^|, it will be increased in 
the outer while loop of the algorithm, until it is sufficiently large. But then, the algorithm 
performs some preliminary computations of the intersection lines and vertices, which cost 
additional run-time. Thus our extrapolation formula is aimed at choosing E slightly greater 
than the exact bound for |;u| we expect. 

Proposition 19. Algorithm^ finds all the hemispheres of radius on which a 2-cell of 
the Bianchi fundamental polyhedron can lie. This algorithm terminates within finite time. 



Proof. 



Directly from the definition of the hemispheres S^^Aj it follows that the radius is 
given by So our algorithm runs through all ^ in question. By construction of 
the Bianchi fundamental polyhedron D, the hemispheres on which a 2-cell of D lies 
must have their centre in the fundamental rectangle Dq. By proposition [131 such 
hemispheres cannot be everywhere below some other hemisphere in the list. 
Now we consider the run-time of the algorithm. There are finitely many ^ £ O the 
norm of which takes a given value. And for a given ^u, there are finitely many X £ O 



such that ^ is in the fundamental rectangle Dq. Therefore, this algorithm consists 
of finite loops and terminates within finite time. 

□ 
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Algorithm 2 Recording the hemispheres of radius jj 

Input: The value A/", and the hemisphere list (empty by default). 
Output: The hemisphere list with some hemispheres of radius jj- added. 

for a running from through M within Z, do 
for b in Z such that \a + buj\ = M, do 
Let fi := a + huj. 

for all the A G O with f in the fundamental rectangle Dq, do 
if the pair (/i, A) is unimodular, then 

Let C be the length of the hemisphere list. 
everywhere_below :— false, j := 1. 
while everywhere_below — false and j ^ £, do 
Let Sp^a be the j'th entry in the hemisphere list; 
if Sfj^^x is everywhere below S'^^q,, then 

everywhere_below :— true, 
end if 

Increase j by 1. 
end while 

if everywhere_below — false, then 

Record S^^x into the hemisphere list, 
end if 
end if 
end for 
end for 
end for 

We recall that the notion "everywhere below" has been made precise in 
definition 1121 : and that the fundamental rectangle Dq has been specified 
in notation |?]. 



Proposition 20. Algorithm^finds the minimal height occurring amongst the non-singular 
vertices ofdB(ai, . . . ,a„). This algorithm erases only such hemispheres from the list, which 
do not change dB(ai, . . . ,a„). It terminates within finite time. 

Proof. 

• The heights of the points in Ti are preserved by the action of the translation 
group Too, so we only need to consider representatives in the fundamental do- 
main Dqc, for this action. Our algorithm computes the entire cell structure of 
dB{ai, . . . , an) D D^o, as described in subsection 16.41 The number of lines to inter- 
sect is smaller than the square of the length of the hemisphere list, and thus finite. 
As a consequence, the minimum of the height has to be taken only on a finite set 
of intersection points, whence the first claim. 

• If a cell of dB{ai, . . . , a„) lies on a hemisphere, then its vertices are lifts of intersec- 
tion points. So we can erase the hemispheres which are strictly below some other 
hemispheres at all the intersection points, without changing dB(ai, . . . ,a„). 
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• Now we consider the run-time. This algorithm consists of loops running through the 
hemisphere list, which has finite length. Within one of these loops, there is a loop 
running through the set of pairs of lines L(^, ^). A (far too large) bound for the 
cardinality of this set is given by the fourth power of the length of the hemisphere 
list. 

The steps performed within these loops are very delimited and easily seen to be of 
finite run- time. 

□ 



Algorithm 3 Computing the minimal proper vertex height 
Input: The hemisphere list {S{ai), . . . , S{an)}- 

Output: The lowest height ^ of a non-singular vertex of dB{ai, . . . , a„). And the hemi- 
sphere list with some hemispheres removed which do not touch the Bianchi fundamental 
polyhedron. 

for all pairs S'^.q,, S^^^x in the hemisphere list which touch one another, do 

compute the line L{^,^) of notation 1101 
end for 

for each hemisphere 5*^ q, in the hemisphere list, do 

for each two lines ^) and L{^,^) referring to ^, do 

Compute the intersection point of L{^, ^) and L{^, ^), if it exists. 

end for 
end for 

Drop the intersection points at which 5^ „ is strictly below some hemisphere in the list. 
Erase the hemispheres from our list, for which no intersection points remain. 
Now the vertices of i?(ai, . . . , a„) Pi D^o are the lifts (specified 

in definition I14D on the appropriate hemispheres of the remaining 
intersection points. 

Pick the lowest value C > for which (z, C) & H is the lift on some hemisphere of a 
remaining intersection point z. 
Return ^. 



7.2. The cell complex and its orbit space. With the method described in subsec- 
tion 16.41 we obtain a cell structure on the boundary of the Bianchi fundamental polyhed- 
ron. The cells in this structure which touch the cusp oo are easily determined: they are four 
2-cells each lying on one of the Euclidean vertical planes bounding the fundamental domain 
Z^oo for Too specified in notation and four 1-cells each lying on one of the intersection 
lines of these planes. The other 2-cells in this structure lie each on one of the hemispheres 
determined with our realization of Swan's algorithm. 

As the Bianchi fundamental polyhedron is a hyperbolic polyhedron up to some missing 
cusps, its boundary cells can be oriented as its facets. Once the cell structure is subdivided 
until the cells are fixed pointwise by their stabilisers, this cell structure with orientation is 
transported onto the whole hyperbolic space by the action of T. 
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7.3. Computing the vertex stabilisers and identifications. Let us state explicitly the 
F-action on the upper-half space model in the form in which we will use it rather than 
in its historical form. 



Lemma 21 (Poincare). If ^ ~\c d)^ GL2(C), the action of j onH is given by 
7 • {z,C) = {z',C'), where 

i_ |det7|C , _ [d - cz) (az - b) - (^ca 

From this operation formula, we establish equations and inequalities on the entries of a 
matrix sending a given point {z,() to another given point {z'X') iii We will use them 
in algorithm U to compute such matrices. For the computation of the vertex stabilisers, we 
have {z,C) = {^',0^ which simplifies the below equations and inequalities as well as the 
pertinent algorithm. First, we fix a basis for O as the elements 1 and 

{\J —m^ m = 1 or 2 mod 4, 

~2 ~'~ 2 V m = 3 mod 4. 

As we have put m 7^ 1 and m 7^ 3, the only units in the ring O are ±1. We will use the 
notations \x\ := min{n € Z | n ^ x} and [xj := max{n E Z | n ^ x} for x G M. 

Lemma 22. Let m = 3 mod 4. Let ^ j € SL2(0) be a matrix sending {z, r) to ((^, p) eH. 



Write c in the basis as j -\- kuj , where j, k £ Z. Then |c|2 ^ |j| ^ ~Tp^ '^^^ 



/m+l _ -2™ ^ . m+1 _ -2™ 

2? y rp 3 m 2j V rp ^ 



m+l m + l m+l m + l 

Proof. From the operation equation ^}\-{z.,r) = {Cp), we deduce \cz — dp +r2|c|2 



and conclude r'^\c\'^ ^ ^, whence the first inequality. We insert |c|2 = ( j — |) + m (|) 
_ j2 _|_ !n±lk'2 — jf^ into it, and obtain 

m + l m + l\ rp J 

We observe that f{k) is a quadratic function in A; G Z C M, taking its values exclusively 
in R. Hence its graph has the shape of a parabola, and the negative values of f(k) appear 
exactly on the interval where k is between its two zeroes, 

k± = ^± 2^, where A = ^ - j2m. 

m+l m+l ' rp 

This implies the third and fourth claimed inequalities. As A; is a real number, A must be 

i+ — 

non-negative in order that f{k) be non-positive. Hence j2 ^ , which gives the second 
claimed inequality. □ 
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Algorithm 4 Computation of the matrices identifying two points in T-L. 



Input: The points {z,r), {C,p) in the interior of T-L, where z, ( G K and r^, G Q. 
Output: The set of matrices G SL2(C'_m), m = 3 mod 4, with nonzero entry c, 

sending the first of the input points to the second one. 

c will run through O with < IcP ^ — . 
Write c in the basis as j + koj , where j,k G Z. 



for j running from 



i+i 

m 

rp 



through 



rp 



do 



k 



limit ■ 



1 



m+1 



±2- 



rp ^ 

m+1 



for k running from L^iTmitJ through [feymitl 
c := j + few; 

if IcP ^ — and c nonzero, then 

Write cz in the basis as R{cz) + W{cz)uj with i?(cz), T^(c2;) G 
d will rim through O with |cz - dp + r^|cp = ^. 
Write d in the basis as q + su, where q,s G Z. 



'limit 



-r2|c|2 



for s running from [siij^j^J through rsi[„jjtl do 

|2 /W(cz) 
\ — m. \ — i — '- 



if A is a rational square, then 



W{cz) 



+ - ± 



g± := R{czj — 2 
Do the following for both q± 
if q± G Z, then 

d := g± + soj; 



q+ and g± = g_ if A 7^ 0. 



if 



\cz 



dp + r^lcp 



and (c, d) unimodular, then 



a:= fd-fcz-cC 
if a is in the ring of integers, then 

6 is determined by the determinant 1 : 
b 



ad—l 



Return 

end if 
end if 
end if 
end if 
end if 
end for 
end if 
end for 
end for 



if b is in the ring of integers, then 
Check that ^)-(z,r) = (C,p). 

d 
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Lemma 23. Under the assumptions of lemma [22l write d in the basis as q + soj, where 
q,s G Z. Write cz in the basis as R{cz) + W{cz)lo, where R{cz),W{cz) G Q. Then 

W{cz) - 2^J^-^ ^ s ^ W{cz) + 2^^-^, and 



^, , Wicz) s r /W(cz) ^ 

R(cz) — + -±\--r^c^-m' ^ ' 

2 2 V P 



Proof. Recall that co = + 2\f~^i sog + sa; = g — | — l^^—m. The operation equation 
yields \cz — dp + r^|cp = ^. From this, we derive 



^2u|2 



^ r-\cy = (cz- (g + sw)) (cz- (g- f - f\/^)) 
= (Re(cz) - g + f )^ + (Im(cz) - f V^)^ 
= Re(c2;)^ + - 9S + X ~ 2Re(cz)g + Re(cz)s + (lm(c2;) - f 

We solve for g, 



2| |2 

r c 



+ (-2Re(c2;) - s) g + (^Re(cz) + |) V (lm(c2:) - i^ypm^ - ^ + 

and find 

g± = Re(cz) + f ± a/A, where A = ^ — r^|cp — (im(cz) — l^y/mf' . 
We express this as 

g± = i?(cz) - ^ + f ± ^/A, where A = ^ - r^jcp - m - f " 

which is the claimed equation. The condition that q must be a rational integer implies 
A ^ 0, which can be rewritten in the claimed inequalities. □ 

We further state a simple inequality in order to prove that algorithm |4] terminates in 
finite time. 



Lemma 24. Let K = Q{^/—m) with m ^ 3. Let c,zGK. Write their product cz in the 
Q'basis {1,^;} for K as R{cz) + W{cz)uj. Then the inequality \ W{cz)\ ^ |c| • \z\ holds. 

Proof. Let x + yu G K with j;, y G Q. Our first step is to show that |y| ^ |x + yu:\. Consider 
the case m= 1 or 2 mod 4. Then 

\x + yuj\ = \/ x"^ + my"^ ^ V^\y\ ^ \y\i 
and we have shown our claim. Else consider the case m = 3 mod 4. Then, 




\x + yu)\ = {x + ojy){x + ojy) = 

and our claim follows for m > 3. Now we have shown that |VF(cz)| ^ \cz\; and we use some 
embedding of K into C to verify the equation \cz\ = \c\ ■ \z\. □ 

Proposition 25. Letm = 3 mod 4. Then algorithm\^ gives all the matrices ^^G SL2(C') 

with c^ 0, sending {z,r) to {(, p) S Ti. It terminates in finite time. 

Proof. 

• The first claim is easily established using the bounds and formulae stated in lemmata 
[22] and [231 
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• Now we consider the run-time. This algorithm consists of three loops the limits 
of which are at most linear expressions in For we use lemma [24] and 

r^|cp ^ J to see this (we get a factor \z\ here, which we can neglect). 

□ 

Finally, it should be said that the scope of computations one can do with geometric 
models for the Bianchi groups does not stop once the integral homology of the full group 
is known. There is further interest in homology with twisted coefficients, congruence sub- 
groups and modular forms (see for instance |§enll| , |§T09| ). Currently, Page [Pagl2 is 
working on optimizing algorithms in order to obtain more cell complexes for Bianchi groups 
and other Kleinian groups. 
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